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INTRODUCTION 


BACKGROUND 

The  Military  Specification  for  Flying  Qualities  of  Piloted  Airplanes,  Ml  L-F-8785B,  (reference 
(a)),  was  developed  largely  from  flight  tests  of  classically  responding  unaugmented  aircraft.  Its 
quantitative  requirements  are  generally  expressed  in  terms  of  modal  approximations  which  can  be 
described  mathematically  by  first  or  second  order  linear  expressions.  Advancements  in  aerody¬ 
namics  and  complicated  control  system  augmentation  schemes,  prevalent  in  modern  aircraft  designs, 
have  resulted  in  responses  which  are  described  by  high  order  functions.  The  high  order  modal 
parameters  may  differ  significantly  from  their  linear  components. 

In  an  attempt  to  utilize  the  existing  requirements  in  analyzing  advanced  aircraft/control 
system  configurations,  the  concept  of  equivalent  systems,  has  been  introduced  (reference  (b)).  A 
digital  frequency  domain  equivalent  system  matching  technique  has  been  developed  by  Hodgkinson, 
et.  al .,  and  applied  to  the  high  order  representations  of  experimental  aircraft  (references  (c)  and  (d)). 
The  approach  used  was  to  approximate  the  high  order  pitch  rate  to  pilot  control  input  transfer 
function  of  the  subject  aircraft  with  a  classical  low  order  transfer  function  describing  the  specifi¬ 
cation  requirements,  augmented  with  a  time  delay.  This  equivalent  time  delay  approximates  the 
phase  lag  introduced  by  the  high  frequency  control  system  components.  Within  the  scope  of  the 
initial  investigations,  it  was  determined  that  the  linear  modal  requirements  of  MIL-F-8785B,  when 
augmented  by  a  requirement  on  time  delay,  are  appropriate  for  specifying  the  handling  qualities 
of  the  advanced  high  order  configurations  of  tomorrow's  airplanes  (reference  (e)).  This  approach 
has  been  incorporated  in  the  latest  revision  to  the  MIL-SPEC,  MIL-F-8785C,  reference  (f),  which 
states: 


"The  contractor  shall  define  equivalent  classical  systems  which  have  responses  most 
closely  matching  those  of  the  actual  aircraft." 

The  parameters  defining  the  resulting  equivalent  system  (frequency,  damping  ratio,  etc.)  rather 
than  any  modes  of  the  high  order  system,  are  to  be  compared  with  the  specification  requirements. 
However,  no  guidance  is  given  as  to  how  the  contractor  shall  proceed  with  his  equivalent  system 
definition  nor  with  what  criteria  its  adequacy  will  be  judged  by  the  procuring  agency. 


The  Naval  Air  Development  Center,  as  part  of  its  effort  in  identifying  flying  qualities  criteria 
for  manned  aircraft,  undertook  the  determination  of  equivalent  system  descriptions  of  current  Navy 
tactical  aircraft.  This  investigation  not  only  encompassed  the  determination  of  the  classical  pitch- 
rate  short-period  model  for  current  fleet  aircraft,  but  also  investigated  the  effects  of  various  high 
order  configuration  components  and  additional  aircraft  response  parameters. 

PURPOSE 

The  purpose  of  this  effort  was  to  investigate  the  utility  of  the  equivalent  systems  approach  to 
defining  the  dynamic  longitudinal  flying  qualities  parameters  of  augmented  aircraft.  This  report 
presents  equivalent  low  order  system  models  for  current  U.S.  Navy  tactical  aircraft  and  compares 
them  with  the  modal  requirements  of  MIL-F-8785C,  reference  (f). 

SCOPE 

The  low  order  equivalent  systems  presented  in  this  report  were  determined  via  frequency 
response  matching  techniques.  The  relative  merits  of  frequency  versus  time  response  matching 
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have  been  addressed  elsewhere  (reference  (g)),  and  will  not  be  further  discussed  in  this  report.  The 
frequency  response  matching  technique  was  chosen  for  this  effort  due  to  its  more  definitive  state  of 
development  and  the  fact  that  it  is  currently  being  investigated  throughout  the  aircraft  industry. 

The  aircraft  included  in  this  analysis  were  the  A-6,  A-7,  F-14,  F-18  and  S-3.  Where  applicable, 
each  aircraft  was  assumed  to  have  its  Stability/Control  Augmentation  System  (SAS/CAS)  ON.  Only 
longitudinal  dynamics  were  analyzed.  The  flight  conditions  investigated  included  both  Power 
Approach  (PA)  and  Cruise  (CR)  configurations  as  presented  in  table  I. 


TABLE  I.  FLIGHT  CONDITIONS 


Aircraft 

Configuration 

Gross 

Weight 

(lb) 

CG 

Position 
(%  MAC) 

Altitude 

(ft) 

Airspeed 

(M/KEAS) 

A-6 

CR 

35905 

23.6 

20000 

0.4 

0.72 

0.88 

PA 

26600 

24.1 

0 

95 

S-3 

CR 

36320 

21.7 

15000 

0,36 

0.71 

PA 

31790 

25.0 

0 

97 

F-14 

CR 

51015 

8.2 

15000 

0.5 

8.2 

0.7 

10.4 

0.83 

11.4 

1.2 

PA(D/o\ 

46380 

9.8 

0 

126 

PA(2)'2' 

9.8 

0 

121 

A-7 

CR 

21890 

30.0 

15000 

0.3 

0.6 

0.9 

F-18 

CR 

29930 

25.0 

10000 

0.5 

PA 

30700 

25.0 

0 

133 

Notes:  (1)  Direct  Lift  Control  (DLC)  ON 
(2)  DLC  OFF 


METHOD 

Frequency  response  matching  techniques  were  utilized  to  determine  low  order  equivalent 
systems  describing  the  complex  aircraft  high  order  representations.  Digital  computer  programs, 
prepared  by  the  McDonnell  Aircraft  Company,  utilized  a  direct  Rosenbrock  search  algorithm  (refer¬ 
ence  c)  to  match  a  Bode  plot  describing  the  high  order  pilot  input  to  aircraft  output  transfer  func¬ 
tion  with  an  equivalent  low  order  system.  Since  this  analysis  is  concerned  with  determining  equiva¬ 
lent  longitudinal  short  period  models,  the  pitch  rate  and  normal  acceleration  responses  to  pilot 
control  inputs  were  analyzed. 
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In  order  to  use  the  matching  routines,  a  description  of  the  frequency  response  of  the  system 
to  be  matched  is  required.  This  may  be  in  the  form  of  either  (1)  transfer  functions  or  (2)  numerical 
phase-gain  data  obtained  at  various  input  frequencies.  Since  only  limited  numerical  response  data  is 
available  for  the  subject  airplanes  (and  it  is  generally  corrupted  with  instrumentation  noise  and  air 
turbulence)  the  transfer  function  input  approach  was  chosen.  Each  aircraft's  transfer  functions 
describing  the  desired  responses  were  obtained  either  directly  from  available  information  (A-7)  or 
computed  via  NADC  transfer  function  programs  from  stability  and  control  derivative  information 
(A-6,  S-3,  F-14,  F-18).  References  (h)  through  (m)  were  used  to  obtain  this  information  as  well  as  a 
description  of  the  respective  control  systems.  With  the  aircraft's  unaugmented  dynamics  thus 
obtained,  the  control  components  present  in  each  aircraft's  control  system  (i.e.,  actuators,  stick 
feed  system,  feedback  loops,  compensation  networks,  etc.)  were  added  to  obtain  the  high  order 
transfer  function  describing  each  aircraft/control  system  combination  and  flight  condition.  Brief 
descriptions  of  the  aircraft  and  their  respective  control  systems  are  presented  in  appendix  A. 


Two  computer  programs,  LONFIT  (reference  n)  and  NAVFIT  (reference  o)  were  utilized  to 
obtain  the  frequency  response  matches  presented  in  this  report.  The  LONFIT  program  utilized  a 
low  order  system  of  the  form: 


K0  (s  +  La)  e 


-r0s 


x  S2  +  2fSD  CO 5q  S  +  coir 


>sp^sp-  ■  wsp 
when  matching  the  pitch  rate  response,  and 
ri 


Knz  e*TnzS 


s2  +  2fSp  coSp  s  + 


sp 


(1) 


(2) 


when  matching  the  normal  acceleration  response.  The  NAVFIT  program  is  a  general  program 
which  allows  the  analyst  to  choose  the  order  of  the  assumed  low  order  system. 


The  match  was  obtained  by  minimizing  the  sum  of  the  squares  of  the  difference  in  both  mag¬ 
nitude  and  phase  of  the  high  and  low  order  system  as  expressed  by: 

M  =  2  ^(Gain|_|Q5  -  Gain^Qg)^  +  0.01745  (Phase^Qg  -  Phase|_Qg)^j  (3) 

where  gain  is  in  decibels  and  phase  in  degrees. 


The  summation  was  performed  at  a  number  of  discrete  intervals  (typically  21 )  over  the  pilot's 
short  period  frequency  range  of  interest  (nominally  .1  to  10  rad/sec). 


In  performing  the  match,  the  analyst  worked  interactively  with  the  program  to  determine 
which  of  the  decision  variables  (Kg,  La,  fSp,  coSp,  Tq)  to  vary  in  the  search  for  the  equivalent  low 
order  system.  In  general,  the  McAl  R  recommended  procedure,  outlined  below,  was  utilized  to 
introduce  decision  variables  into  the  search. 


a.  Initially,  attempt  to  match  the  high  order  system  response  in  the  frequency  range  of  0.1 
to  10  rad/sec  by  varying  gain  or  zeta  or  omega  or  all  three,  with  La  fixed  at  the  estimated 
airplane  value,  and  with  no  time  delay. 

b.  If,  after  obtaining  a  plot  of  the  high  order  system  gain  characteristics,  it  is  evident  that 
the  0.1  through  10  rad/sec  frequency  range  does  not  include  the  peak  in  gain,  expand 
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the  frequency  range  to  include  the  peak.  If  two  separate  peaks  are  evident  in  the  0.1 
through  10  rad/sec  frequency  range,  modify  the  range  to  span  only  the  assumed  short 
period  peak. 

c.  If  a  good  match  (a  payoff  function  less  than  20  is  considered  a  good  match)  cannot  be 
obtained  after  following  steps  a  and  b,  include  time  delay  in  the  search.  (Inclusion  of 
the  time  delay  is  usually  effective  when  a  mismatch  in  phase  at  higher  frequencies  is 
evident.) 

d.  If  after  following  steps  a,  b  and  c,  a  low  payoff  has  not  resulted,  delete  the  time  delay 
and  include  La  in  the  search.  (Inclusion  of  !_a  is  usually  effective  when  a  mismatch 
in  low  and  intermediate  frequencies  is  evident.) 

e.  If  after  following  steps  a,  b,  c,  and  d,  a  good  match  has  not  been  obtained,  include  both 
La  and  time  delay  in  the  search. 

f.  If  after  the  above  procedure,  an  acceptable  payoff  function  cannot  be  obtained,  the 
configuration  dynamics  must  be  considered  to  not  be  represented  by  the  assumed  low 
order  system  and  an  alternative  evaluation  method  is  necessary. 

In  performing  these  matches,  the  "goodness  of  fit"  was  determined  quantitatively  by  the  value 
of  the  mismatch  parameter  and  qualitatively  from  Bode  plots  and  time  history  responses  to  unit 
impulse  and  step  control  inputs.  Frequency  and  time  history  comparisons  for  representative  equi¬ 
valent  systems  examples  are  presented  in  appendix  B. 

The  resulting  modal  parameters  representing  the  low  order  systems,  as  well  as  the  dominant 
roots  of  the  high  order  system,  were  compared  against  the  requirements  of  MIL-F-8785C.  The 
dominant  roots  of  the  high  order  system  were  defined  as  that  oscillatory  pair  present  in  the  fre¬ 
quency  range  from  0.5  to  10  rad/sec. 


RESULTS  AND  DISCUSSION 


GENERAL 

Equivalent  system  models  were  initially  obtained  for  the  pitch  rate  to  cockpit  control  position 
transfer  functions  for  the  A-6,  S-3  and  F-14  airplanes.  In  these  cases,  the  matching  procedure  was 
straightforward,  resulting  in  excellent  frequency  and  time  history  response  matches.  Subsequently, 
control  stick  dynamics  were  included  in  the  F-14  airplane's  high  order  system  model  and  equivalent 
pitch  rate  to  control  force  transfer  functions  were  generated.  Control  force  was  also  used  as  the 
input  command  for  the  A-7  and  F-18  airplanes,  although  in  these  cases,  control  stick  dynamics  were 
not  modeled. 

The  inclusion  of  feel  system  dynamics  compounded  the  matching  problem  by  introducing 
additional  roots  within  the  short  period  frequency  range  of  interest.  As  a  result,  these  equivalent 
systems  resulted  in  either  large  mismatch  parameters  or  large  numerator  root  (La)  values  which  are 
inconsistent  with  the  airplane's  lift  curve  slope. 

Two  methods  of  restricting  the  variation  in  La  were  analyzed.  First,  the  normal  acceleration 
response  was  included  into  the  classical  short  period  approximation  matching  procedure.  Matching 
the  high  order  normal  acceleration  at  the  instantaneous  center  of  rotation  provided  a  short  period 
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characteristic  equation  consistent  with  that  obtained  from  the  La  fixed  pitch  rate  analysis.  In 
addition,  simultaneously  matching  pitch  rate  and  normal  acceleration  significantly  reduced  the 
variation  in  equivalent  La.  Secondly,  since  the  majority  of  the  mismatch  evident  with  La  fixed 
occurred  in  the  vicinity  of  the  additional  feel  system  roots,  the  equivalent  short  period  definition 
was  modified  to  account  for  this  additional  root.  Matching  pitch  rate  and  normal  acceleration  with 
first  over  third  and  zero  over  third  order  transfer  functions,  respectively,  resulted  in  low  mismatch 
values  and  excellent  time  history  comparisons. 

Frequency  and  time  history  comparisons  for  each  of  the  identified  equivalent  systems  models 
are  presented  in  appendix  B. 

Comparison  of  the  equivalent  time  delay  and  damping  ratio  with  the  requirements  of  Ml  L-F- 
8785C  were  straightforward.  This  may  be  attributed  to  the  fact  that  the  requirements  are  imposed 
singularly  on  the  parameters  of  interest  and  are  not  correlated  with  any  other  parameters.  The  short 
period  frequency  requirement,  however,  correlates  frequency  (tosp)  with  acceleration  sensitivity 
(n/a)  via  boundaries  of  constant  control  anticipation  parameter  (a?2sp/n/a).  It  is  a  simple  matter  to 
plot  equivalent  short  period  frequency  as  a  function  of  equivalent  acceleration  sensitivity.  However, 
the  resulting  classical  equivalent  short  period  control  anticipation  parameter  does  not  correlate 
with  the  high  order  system's  control  anticipation  parameter  (#rnax/nzss)*  To  alleviate  this  problem 

a  control  anticipation  parameter  was  developed  which  accounts  for  the  attenuation  affects  of  the 
high  order  control  system  (reference  (p)).  The  use  of  this  parameter  results  in  a  consistent  descrip¬ 
tion  of  both  the  high  and  low  order  systems  which  can  be  correlated  to  pilot  opinion  boundaries. 

CLASSICAL  SHORT  PERIOD  EQUIVALENT  SYSTEMS 

The  short  period  mode  of  motion  of  interest  is  primarily  characterized  by  a  rapid  rotation  in 
aircraft  attitude.  Therefore,  the  major  emphasis  in  longitudinal  equivalent  systems  matching  has 
been  centered  on  the  pitch  response,  the  parameters  of  which  are  necessary  for  comparison  with 
the  requirements  of  Ml  L-F-8785C.  The  same  approach  has  been  taken  in  this  analysis. 

Initial  matches  were  obtained  for  the  pitch  rate  to  pilot  control  position  input  transfer  func¬ 
tions  for  the  A-6,  S-3  and  F-14  airplanes.  In  general,  these  configurations  were  characterized  by 
the  bare  airframe  augmented  by  control  surface  actuators  and  pitch  rate  feedbacks.  From  these 
configurations,  a  basic  understanding  of  the  frequency  response  matching  process  was  obtained. 
Subsequently,  the  effect  of  compounding  dynamic  components  such  as  control  stick  dynamics, 
normal  acceleration  feedbacks  and  system  prefilters  were  investigated  by  matching  pitch  rate  to 
pilot  force  inputs  for  the  F-14,  A-7  and  F-18  aircraft. 

A-6  AIRPLANE 

The  A-6  airplane  utilizes  washed  out  pitch  rate,  as  described  in  appendix  A,  to  augment  the 
airplane's  basic  stability  characteristics.  The  resulting  pitch  rate  response  to  pilot  longitudinal 
control  position  inputs  can  be  represented  by  a  fourth  order  numerator  over  sixth  order  denom¬ 
inator  transfer  function.  The  pitch  rate  feedback  component  adds  a  numerator  and  denominator 
root  in  the  vicinity  of  0.5  rad/sec  while  the  actuator  adds  a  single  denominator  root  at  approxi¬ 
mately  30  rad/sec.  The  proximity  of  the  added  feedback  roots  results  in  their  effectively  cancelling 
one  another  while  the  actuator  root  is  quite  far  removed  from  the  closed  loop  short  period  root,  it 
is  to  be  expected  therefore,  that  the  equivalent  roots  identified  by  the  frequency  matching  pro¬ 
cedure  will  quite  closely  match  the  "dominant  roots"  of  the  high  order  system  augmented  with  a 
time  delay. 
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The  McAir  recommended  matching  procedure  was  utilized  in  the  NAVFiT  program  to  deter¬ 
mine  a  first  order  numerator  over  second  order  denominator  equivalent  system  for  each  of  the  flight 
conditions  analyzed.  The  complete  procedure  was  followed,  even  though  acceptable  matches  were 
obtained  by  varying  only  gain,  damping  ratio  and  frequency,  in  order  to  obtain  an  understanding  of 
the  effects  of  introducing  different  decision  variables  into  the  search.  The  A-6  high  and  low  order 
equivalent  systems  are  summarized  in  table  II. 

With  the  phugoid  roots  included  in  the  high  order  model,  the  frequency  range  over  which  the 
match  was  conducted  had  to  be  modified  to  exclude  their  contributions.  The  effect  of  varying 
the  match  frequency  range  to  exclude  phugoid  contributions  is  presented  in  figure  1.  With  the 
frequency  matching  range  including  the  tail  of  the  phugoid  resonant  peak,  the  equivalent  short 
period  damping  and  frequency  are  different  from  that  of  the  high  order  system  dominant  root 
pair  and  a  high  mismatch  results. 

Modifying  the  frequency  range  to  exclude  the  phugoid  contributions  yields  results  which 
correspond  with  those  obtained  from  removing  the  phugoid  roots  from  the  high  order  system  prior 
to  performing  the  match.  In  the  cases  in  which  the  phugoid  contributions  were  not  excluded,  the 
search  routine  is  attempting  to  match  the  gain  and  phase  characteristics  of  two  resonant  peaks 
with  only  one  set  of  numerator  and  denominator  breakpoints.  As  a  result,  a  compromise  is 
achieved  which  provides  the  best  average  mismatch  across  the  entire  frequency  range,  but  only 
approximates  the  desired  short  period  characteristics. 

In  the  majority  of  the  cruise  configurations  analyzed,  the  phugoid  and  short  period  roots 
were  widely  enough  separated  that  modifying  the  match  frequency  range  produced  acceptable 
results.  However,  in  the  power  approach  cases,  the  phugoid  and  short  period  roots  were  much 
more  closely  coupled.  In  these  instances,  increasing  the  low  frequency  boundary  to  exclude  phugoid 
effects  also  cut  off  some  of  the  lower  frequency  short  period  information.  For  the  A-6  power 
approach  case,  the  lower  frequency  bound  was  raised  to  0.5  rad/sec  to  exclude  phugoid  contribu¬ 
tions.  However,  the  aircraft's  numerator  term,1/Tg2#  is  at  0.461  rad/sec.  Therefore,  using  0.5 
rad/sec  as  the  lower  frequency  bound  cuts  off  a  portion  of  the  short  period,  as  well  as  the  phugoid 
response.  Although  acceptable  matches  were  attainable  by  using  a  lower  frequency  bound  of  0.5 
rad/sec,  the  resulting  short  period  parameters  will  later  be  shown  to  have  less  than  satisfactory 
correlation  with  the  high  order  modal  parameters.  In  instances  such  as  this,  where  the  phugoid  and 
short  period  roots  were  closely  coupled,  the  phugoid  roots  were  removed  from  the  high  order 
model  prior  to  performing  the  equivalent  system  match. 

With  the  frequency  range  modified  to  exclude  the  phugoid  contributions  to  the  gain  and  phase 
plots,  the  major  improvement  in  system  match  was  obtained  via  the  time  delay  parameter.  The 
time  delay  serves  to  modify  the  phase  characteristics  of  the  low  order  system  to  account  for  the 
phase  characteristics  attributable  to  the  high  frequency  terms  in  the  high  order  system.  It  has 
unity  gain  and  therefore  affects  only  the  Bode  phase  characteristics.  Inclusion  of  the  time  delay 
parameter  had  little  effect  on  the  identified  frequency  and  damping  ratios  but  greatly  reduced  the 
mismatch  parameter  by  improving  the  phase  match  at  the  higher  frequencies,  as  shown  in  figure  2. 
The  equivalent  pitch  rate  time  history  is  also  improved,  although  there  is  no  response  from  t  =  0  to 
t  =  r  and  the  initial  value  of  pitch  acceleration  is  much  higher  than  that  of  the  high  order  system 
in  order  to  match  the  slope  of  the  pitch  rate  response  following  the  delay. 

Freeing  l_a  in  the  search  process  also  improved  the  quality  of  the  match,  although  generally 
not  as  much  as  the  time  delay  parameter  did.  As  can  be  seen  from  figure  3,  freeing  La  primarily 
improved  the  match  at  frequencies  below  the  short  period  frequency  and  only  slightly  improved 
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TABLE  II.  A-6  AIRPLANE  EQUIVALENT  TRANSFER  FUNCTIONS 
Pitch  Rate  Response  to  Cockpit  Control  Position  Inputs 
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Note:  (1)  Phugoid  contributions  removed  from  high  order  system 
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the  match  at  the  higher  frequencies.  Freeing  La,  with  zero  time  delay,  resulted  in  higher  values 
of  La  than  freeing  both  La  and  time  delay  due  to  the  attempt  at  matching  the  high  frequency 
phase  characteristics  with  fewer  parameters. 

Based  on  these  results,  phugoid  contributions  were  excluded  from  the  matching  process,  time 
delay  was  included  in  all  subsequent  matches,  and  La  was  allowed  to  vary,  where  necessary,  to 
improve  the  mismatch  parameter. 

S-3  AIRPLANE 

There  are  no  stability  augmentation  components  included  in  the  longitudinal  control  system 
of  the  S-3  airplane.  Its  high  order  representation  consists  of  unaugmented  phugoid  and  short 
period  dynamics  controlled  by  an  elevator  actuator  with  a  0.0294  sec.  time  constant. 

The  NAVFIT  frequency  matching  routine  readily  identifies  the  S-3's  short  period  dynamics 
with  equivalent  delays  of  0.027—0.029  sec.  as  presented  in  table  III.  The  power  approach  charac¬ 
teristics  again  necessitate  the  removal  of  phugoid  dynamics  from  the  high  order  representation,  due 
to  the  close  coupling  between  short  period  and  phugoid  roots. 

F-14  AIRPLANE 

Pitch  Rate  Response  Matching  —  The  F-14  airplane's  control  system  includes  a  washed  out 
pitch  rate  feedback  signal  that  is  passed  through  a  shaping  network  to  provide  the  desired  response 
to  command  inputs.  The  resulting  F-14  high  order  pitch  rate  to  pilot  control  position  transfer  func¬ 
tion  can  be  represented  by  a  sixth  order  numerator  over  eighth  order  denominator. 

With  L  *  fixed  at  the  value  obtained  from  the  F-14  airplane's  aerodynamic  characteristics, 
acceptable  low  order  equivalent  systems  (mismatch  ^  10)  are  readily  obtained,  as  presented  in 
table  IV.  This  may  be  attributed  to  the  observation  that  although  there  are  a  large  number  of  roots 
present  in  the  frequency  range  of  interest,  each  closed  loop  numerator  root  is  accompanied  by  a 
similar  denominator  root  and  therefore  only  effects  the  response  over  a  limited  range  of  fre¬ 
quencies. 

It  can  also  be  seen  that  freeing  La'  results  in  a  significant  improvement  in  the  mismatch  param¬ 
eter  to  values  less  than  1.  This  reduction  in  mismatch  is  accompanied  by  increases  in  La  and  re¬ 
ductions  in  time  delay,  damping  ratio  and  gain.  The  frequency  increases  for  subsonic  flight  condi¬ 
tions  and  decreases  for  the  one  supersonic  case. 

The  feel-system  in  the  F-14  airplane  includes  a  sprashpot  to  damp  out  control  stick  oscillations 
and  a  feel  spring  to  provide  pilot  force  feel  of  commanded  inputs.  This  feel  system  may  be  repre¬ 
sented  by  a  first  over  third  order  transfer  function  as  shown  in  appendix  A.  The  addition  of  this 
feel  system  to  the  augmented  F-14  airplane's  pitch  dynamics  results  in  a  seventh  order  numerator 
over  eleventh  order  denominator  pitch  rate  transfer  function. 

The  equivalent  systems  obtained  for  the  stick  force  command  inputs  are  included  in  table  IV. 
Matching  the  stick  force  command  high  order  transfer  functions  with  first  over  second  order  equiva¬ 
lent  systems  resulted  in  significantly  different  modal  parameters  from  those  obtained  in  the  control 
position  analysis.  With  L  fixed,  the  equivalent  frequency  is  lower  and  the  time  delay  and  mis¬ 
match  parameters  are  considerably  higher. 
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Note:  (1)  Phugoid  contributions  removed  from  high  order  system 
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Note:  (1)  Phugoid  components  removed  from  high  order  system  prior  to  performing  match 
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The  pitch  response  resulting  from  step  control  and  force  inputs  are  compared  in  figure  4, 
where  the  magnitude  of  the  inputs  has  been  normalized  to  yield  the  same  steady  state  pitch  rate. 
The  addition  of  the  control  stick  feel  system  results  in  a  more  sluggish  initial  pitch  response  than 
that  obtained  for  control  position  inputs  (i.e.— the  frequency  is  lower  and  the  response  is  delayed— 
as  indicated  by  comparing  the  force  and  position  equivalent  system  parameters). 


CONTROL  POSITION  INPUT 


CONTROL  FORCE  INPUT 


Figure  4  —  F-14  Airplane  Response  to  Control  Force  and 
Position  Commands  —  .5M/1 5,000  ft  Altitude 

Freeing  La  reduces  the  mismatch  parameter  to  values  similar  to  those  obtained  in  the  control 
position  analysis.  This  is  accomplished  at  the  expense  of  La;  however,  which  is  now  seen  to  "gallop" 
to  extremely  large  values  in  an  attempt  to  match  the  gain  characteristics  of  the  high  order  system. 
Where  freeing  La  in  the  control  position  cases  resulted  in  increases  in  La  of  30  to  75  percent,  it  now 
increases  by  75  to  500  percent.  The  resulting  frequency  and  time  history  responses  show  very  good 
agreement  between  the  high  and  low  order  systems  (Appendix  B,  figure  B~6).  However  it  is  not 
conceivable  that  the  control  system  implementation  could  alter  La  to  the  extent  that  the  equivalent 
system  analysis  would  indicate. 

Normal  Acceleration  Response  Matching  —  The  equivalent  systems  obtained  from  the  pitch 
rate  transfer  function  often  resulted  in  conditions  of  galloping  La.  This  was  most  evident  for  those 
cases  in  which  the  control  system  introduced  additional  roots  in  the  pilot's  frequency  range  of 
interest.  The  LON  FIT  program  includes  the  capability  of  matching  aircraft  normal  acceleration 
via  the  equivalent  system  presented  in  equation  (2). 
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The  normal  acceleration  thus  defined  may  be  viewed  as  a  measure  of  the  aircraft's  path 
response  represented  as: 


5e„ 

e 

s 

nz* 

5e 

y+i/t7 

where  the  path  response  lag  parameter,  1/T^,  may  be  shown  to  be  approximately  equal  to  the 
numerator  term  in  the  d/8e  transfer  function.  Therefore  matching  the  nz/5e  transfer  function  should 
provide  information  concerning  the  short  period  characteristic  equation  without  having  to  be  con¬ 
cerned  about  matching  the  numerator  roots. 

Equation  (2)  is,  however,  only  valid  for  normal  accelerations  measured  at  the  airplane's  center 
of  rotation.  At  that  location,  the  initial  vertical  acceleration  of  the  tail  produced  by  a  step  elevator 
input  is  balanced  by  that  due  to  the  aircraft's  pitching  acceleration.  At  any  other  location,  this  con¬ 
dition  does  not  hold  and  two  nonminimum  phase  numerator  roots  are  introduced  in  the  normal 
acceleration  transfer  function: 


Matching  normal  acceleration  at  the  center  of  gravity,  or  the  pilot's  location,  with  a  zero  over 
second  (0/2nd)  order  transfer  function  will  only  approximate  the  high  order  system  representation. 
For  example  consider  the  F-14  airplane  response  to  control  position  inputs  at  0.5M,  15,000  ft. 
presented  in  figure  5. 

The  zero  over  second  order  equivalent  system  only  approximates  the  high  order  system  normal 
acceleration  at  the  center  of  gravity.  Although  the  mismatch  parameter  is  reasonable,  the  short 
period  frequency,  damping  ratio  and  equivalent  time  delay  are  significantly  different  from  those 
previously  obtained  from  the  pitch  rate  transfer  function.  In  addition,  the  equivalent  time  response 
does  not  match  the  initial  reversal  in  normal  acceleration  arising  from  the  nonminimum  phase 
numerator  roots.  The  normal  acceleration  at  the  center  of  rotation  is,  however,  matched  at  least 
as  well  as  the  pitch  rate  expression  and  results  in  identical  frequency  and  damping  ratio  and  similar 
time  delays. 

Simultaneously  matching  pitch  rate  and  normal  acceleration  at  the  center  of  rotation,  with  the 
denominators  constrained  to  be  identical,  results  in  a  slight  improvement  in  the  pitch  rate  match  at 
the  expense  of  the  normal  acceleration  match.  This  technique  was  applied  to  the  F  14  cases  in 
which  increased  by  approximately  100  percent  or  more  when  freed  in  the  search  routine.  The 
high  order  representations  and  resulting  equivalent  systems  are  presented  in  table  V. 

The  equivalent  systems  obtained  for  the  normal  acceleration  response  at  the  center  of  rotation 
are,  in  all  cases,  consistent  with  those  obtained  for  the  pitch  rate  response  with  La  fixed.  Simul¬ 
taneously  matching  the  pitch  rate,  with  La  free,  and  normal  acceleration  transfer  functions  has  the 
effect  of  restricting  the  variation  in  evident  when  matching  pitch  rate  alone. 
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Figure  5  —  Comparison  of  Normal  Acceleration  Equivalent  Responses 
F-14  Airplane  —  .5M/15,000  ft-Cockpit  Control  Position  Input 
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Notes:  (1)  15,000  ft  altitude 

(2)  Sea  Level 

(3)  Phugoid  component  removed  from  high  order  system  prior  to  performing  match 
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A-7  Airplane  —  The  A-7  airplane  contains  both  a  stability  augmentation  system  which  utilizes 
pitch  rate  and  normal  acceleration  feedbacks,  and  a  control  augmentation  system  (CAS)  which  elec¬ 
trically  feeds  force  commands  forward  to  be  summed  with  the  mechanical  command  inputs.  Both 
the  normal  acceleration  and  CAS  signals  are  passed  through  a  prefilter  to  eliminate  high  frequency 
inputs  as  outlined  in  appendix  A.  The  high  order  pitch  rate  to  pilot  force  command  transfer  func¬ 
tion  can  be  represented  by  a  fourth  order  numerator  over  sixth  order  denominator.  The  A-7 
equivalent  system  results  are  summarized  in  table  VI. 

Excellent  equivalent  system  matches  are  obtained  for  the  0.6  and  0.9M  cases  as  evidenced  by 
the  very  low  mismatch  parameter.  The  0.3M  case  on  the  other  hand  exhibits  a  large  mismatch 
parameter  attributable  to  the  large  separation  in  closed  loop  numerator  and  denominator  roots 
introduced  by  the  control  system.  At  the  higher  speeds,  these  roots  have  migrated  toward  one 
another  and  have  little  influence  on  the  resulting  frequency  response. 

The  response  at  0.3  Mach  is  characterized  by  relatively  high  mismatch  and  galloping  La; 
similar  to  the  F-14  control  force  command  responses.  Therefore,  the  normal  acceleration  response 
at  the  center  of  rotation,  was  included  in  the  matching  process,  the  results  of  which  are  included  in 
table  VI.  As  with  the  F-14  cases,  the  normal  acceleration  match  parameters  are  consistent  with 
those  of  the  L^-FIXED  pitch  rate  analysis  and  simultaneous  pitch  rate  and  normal  acceleration 
matching  restrict  the  variation  in  La. 

F  18  Airplane  —  The  F~18  airplane  possesses  a  highly  complex  digital  flight  control  system.  It 
incorporates  numerous  compensated  feedbacks,  stick  shaping,  lead-lag  filters,  etc.,  and  has  separate 
control  law  configurations  for  both  cruise  and  power  approach  flight  conditions  as  outlined  in 
appendix  A.  The  pitch  rate  response  to  control  force  inputs  is  described  by  a  14th  order  transfer 
function  in  cruise  configuration  and  by  an  1 1th  order  transfer  function  in  power  approach. 
Although  the  transfer  functions  are  of  relatively  high  order,  the  individual  numerator  and  denomi¬ 
nator  roots  in  the  short  period  frequency  range  are  of  similar  magnitude  and  have  only  local  influ¬ 
ences  on  the  total  frequency  response.  As  a  result,  very  good  match  statistics  are  obtained  for  the 
low  order  systems  as  presented  in  table  VII.  There  are,  however,  two  highly  damped  oscillatory 
denominator  root  pairs  in  the  cruise  configuration  —  one  arising  from  the  aircraft's  short  period 
and  one  from  the  feedback  network.  It  is  difficult  to  determine,  from  the  high  order  representa¬ 
tion,  how  the  combination  of  these  roots  will  affect  the  aircraft's  response  and  the  pilot's  opinion 
of  it. 


Matching  the  F-18's  cruise  configuration  high  order  pitch  rate  response  with  the  classical  short 
period  approximate  system  results  in  excellent  agreement  between  the  two  systems.  The  resulting 
equivalent  frequency  (cu$p e  =  3.05  rad/sec)  lies,  on  a  Bode  plot,  midway  between  the  two  oscilla¬ 
tory  pairs  evident  in  the  high  order  system. 

The  power  approach  configuration  results  are  similar  to  those  obtained  for  the  other  aircraft. 
Freeing  La  to  improve  the  match  statistics  results  in  a  good  time  history  comparison  but  with  a 
large  value  for  La.  The  variation  in  La  may  be  restricted  by  simultaneously  matching  pitch  rate 
and  normal  acceleration. 
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TABLE  VI 

A-7  AIRPLANE  EQUIVALENT  TRANSFER  FUNCTIONS 
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MODIFIED  EQUIVALENT  SYSTEM  DEFINITION 

The  equivalent  systems  obtained  for  the  F-14  airplane's  response  to  force  commands  were 
characterized  by  high  mismatch  values.  Referring  to  the  Bode  plots  presented  in  appendix  B,  this 
mismatch  is  most  evident  in  the  frequency  range  from  1  to  5  rad/sec.  The  cause  of  this  mismatch 
may  be  traced  to  the  single  feel  system  denominator  root  at  3.37  rad/sec.  Since  there  is  no  numer¬ 
ator  root  which  approximates  this  term,  it  affects  a  large  portion  of  the  frequency  response  curves. 
Attempting  to  match  such  a  high  order  system  with  the  classical  short  period  approximation  re¬ 
stricts  the  quality  of  the  match.  There  are  not  enough  degrees  of  freedom  for  the  search  routine  to 
account  for  the  significant  breakpoints  in  the  frequency  range  of  interest.  Increasing  the  order  of 
the  equivalent  system  to  account  for  such  breakpoints  should  conceivably  improve  the  match 
statistics. 

The  effect  of  accounting  for  significant  components  in  the  high  order  transfer  function  was 
investigated  by  matching  the  F-14  force  command  configurations  with  an  equivalent  transfer 
function  of  the  form: 

e(s)  K0(s+La)e'TS  (5) 

F (s)  (s2  +  2?sp  wsps  +  o;2sp>  (s  +  Pe) 

where Pe  is  the  equivalent  denominate r  root  added  to  account  for  the  uncompensated  pole  evident 
in  the  high  order  transfer  function.  The  results  of  this  matching  process  are  presented  in  table  VIII. 
Inclusion  of  the  additional  root  in  the  equivalent  system  denominator  greatly  improved  the  mis¬ 
match  parameter  of  the  low  order  system;  to  the  point  that  it  is  not  necessary  to  even  consider 
freeing  La.  The  short  period  frequency  is  now  greater  than  that  obtained  from  the  high  order  sys¬ 
tem's  dominant  roots.  The  time  delay  has  been  reduced  to  a  value  similar  to  that  obtained  from  the 
control  position  analysis  indicating  that  a  large  portion  of  the  time  delay  associated  with  the  first 
over  second  order  equivalent  models  may  be  attributable  to  this  additional  feel  system  root. 

The  first  over  second  and  first  over  third  order  equivalent  system  responses  are  compared  to 
the  high  order  system  in  figure  6.  The  additional  denominator  root  improves  the  transient  response 
by  eliminating  the  large  time  delay  thereby  improving  the  initial  pitch  acceleration  characteristics. 
However,  the  assumed  low  order  model  is  now  inconsistent  with  that  from  which  the  specification 
requirements  were  generated.  It  would  therefore  seem  necessary  to  determine  new  specification 
requirements  which  would  account  for  the  additional  denominator  root.  This  condition  would, 
however,  not  be  in  concert  with  the  initial  desire  to  utilize  existing  specification  formats  and 
requirements. 

COMPARISON  WITH  SPECIFICATION  REQUIREMENTS 

The  equivalent  system  parameters  obtained  for  each  of  the  aircraft  were  compared  against  the 
requirements  of  Ml  L-F-8785C  as  summarized  in  table  IX. 

The  damping  ratio  and  time  delay  resulting  from  the  equivalent  system  response  can  be  com¬ 
pared  directly  against  the  requirements  of  Ml  L-F-8785C. 
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TABLE  VIII 

MODIFIED  EQUIVALENT  SYSTEM  MODEL  TRANSFER  FUNCTIONS 
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Notes:  (1)  La  FIXED 

(2)  15,000  ft  Altitude 

(3)  Sea  Level 

(4)  Phugoid  component  removed  from  HOS  prior  to  performing  match 
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Figure  6  —  Equivalent  System  Model  Comparisons  —  F-14  Airplane  —  .5M/1 5,000  ft  Altitude 

Control  Force  Command  Inputs 
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TABLE  IX 

Ml  L-F-8785C  DYNAMIC  LONGITUDINAL  STABI LITY  REQUIREMENTS 


Para.  3.2.2. 1.1 

C0sp2/n/a 

Para.  3  2.2.1, 2 
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Para.  3  5.3 
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— 
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- 
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- 
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Notes:  ( 1 )  cosp  >  1 .0  rad/sec  for  n/a  <  3.5 

(2)  coSp  >  0.6  rad/sec  for  n/a  <  2.25 

(3)  cosp  >  0.85  rad/sec  for  2.5  <  n/a  <  4.5 

(4)  toSp  >  0.6  rad/sec  for  1 .6  <  n/a  <  3.75 


In  order  to  compare  the  short  period  frequency  requirements  against  the  specification,  it  is 
necessary  to  determine  the  equivalent  acceleration  sensitivity  (n/ae),  It  can  be  shown,  under  the 
assumption  of  constant  speed  equations,  that  for  the  steady  state  response  to  elevator  inputs: 

n  An/Se  _  V_  1  ^  _V_  . 

a  A a/8e  9  T  ^  9  a‘ 

Therefore,  an  equivalent  n/a  can  be  obtained  by  multiplying  the  equivalent  La  resulting  from  the 
matching  process  by  V/g. 

The  parameters  resulting  from  the  equivalent  system  matches,  along  with  the  higher  order  sys¬ 
tem  "dominant  root"  parameters,  are  compared  against  the  specification  requirements  in  figures  7 
and  8  for  configurations  CR  and  PA,  respectively.  A  number  of  observations  can  be  made  from 
reference  to  this  data: 

1.  The  frequency  and  damping  ratio  of  the  A-6  and  S-3  "dominant"  high  order  system  roots 
are  consistent  with  the  equivalent  system  parameters.  This  may  be  attributed  to  the  lack 
of  control  compensation  roots  in  their  high  order  systems. 

2.  The  time  delay  parameter  reflects  level  1  flying  qualities  for  all  but  the  F-14  response  to 
force  command  inputs  in  which  level  2  time  delays  result. 

3.  Damping  ratio  is  in  the  level  1  region  for  all  conditions  analyzed. 

4.  The  A-6,  S-3  and  F-18  airplanes  and  the  F-14's  response  to  control  position  inputs  all 
yield  level  1  frequency  characteristics.  The  frequency  value  obtained  from  the  dominant 
root,  L^-fixed  and  La-free  analyses  may,  however,  be  significantly  different, 

5.  The  A-7  and  the  F-14  responses  to  control  force  inputs  not  only  result  in  different  values 
of  frequency  for  the  differing  analysis  techniques,  but  the  identified  parameters  tend  to 
cross  over  flying  qualities  boundaries.  For  example,  the  dominant  F-14  root  at  0.5M/ 
15,000  feet  indicates  level  1  frequency  characteristics,  the  0-La  fixed  and  simultaneous 

6  and  nz  equivalent  systems  indicate  level  2-3  characteristics,  and  the  0-La  free  case  pro¬ 
vides  an  indication  that  the  handling  qualities  are,  in  this  instance,  worse  than  level  3. 
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The  basic  premise  under  which  equivalent  systems  have  been  developed  hypothesizes  that:  an 
equivalent  system  which  has  similar  dynamics  (i.e.  frequency  and  time  history  characteristics)  as  a 
higher  order  system  will  be  similarly  evaluated  by  the  pilot  during  maneuvering  tasks.  However,  the 
final  observation  made  above  indicates  differences  in  flying  qualities  levels  and  therefore  pilot  ac¬ 
ceptability  for  various  low  order  descriptions  of  a  particular  high  order  model.  The  frequency  and 
time  history  comparisons  (Appendix  B)  and  their  respective  mismatch  parameters  are  within  reason¬ 
able  bounds  as  defined  in  previous  investigations.  (It  was  determined,  from  the  flight  tests  of  ref¬ 
ence  q  that  configurations  with  mismatch  parameters  as  high  as  200  were  not  noticeable  to  pilots 
when  evaluating  their  equivalence.)  Therefore,  it  becomes  necessary  to  further  investigate  the  short 
period  frequency  specification  methodology  of  MIL-F-8785C. 

CONTROL  ANTICIPATION  PARAMETER  CORRELATION 


The  longitudinal  short  period  frequency  requirements  of  Ml L  F-8785C  are  not  only  presented 
as  a  function  of  n/a  but  are  also  correlated  with  pilot  opinion  (i.e.  flying  qualities  levels)  by  the  con¬ 
trol  anticipation  parameter  (CAP).  Bihrle  (reference  r)  defined  CAP  as  relating  the  two  responses 
of  primary  interest  to  the  pilot  during  a  pullup  —  the  initial  pitch  acceleration  and  the  steady  state 
normal  acceleration.  He  further  showed  that  for  the  constant  speed  short  period  approximation: 

CAP  =  -  iVL  (7) 

n7  n/a 

^$s 

This  parameter  provides  boundaries  on  pilot  opinion  of  acceptable  short  period  frequency  charac¬ 
teristics  when  presented  in  the  format  of  Ml L-F-8785C. 


The  equivalent  system  model  is  the  same  as  the  short  period  approximation.  Equation  (7)  is, 
therefore,  identically  correct  for  the  equivalent  system  model.  It  is  possible  to  plot  coSpe  vs  n /aQ 
and  correlate  it  with  an  equivalent  control  anticipation  parameter,  CAPe.  However,  in  the  case  of 
higher  order  systems,  it  has  been  shown  by  Di Franco,  reference  (s),  that  the  initial  pitch  accelera¬ 
tion  is  zero  and  builds  to  a  maximum  at  some  time  greater  than  t  =  0+.  In  this  instance,  the  pitch 
acceleration  of  interest  in  defining  the  control  anticipation  parameter  is  the  maximum  pitch  accel¬ 
eration  and  not  the  initial  pitch  acceleration.  The  control  anticipation  parameter  is  then  defined  as: 


CAP' 


h 

°max 


0nd 


(8) 


where  0ncj  is  a  non-dimensional  pitch  acceleration.  It  is  the  ratio  of  pitch  acceleration,  including 
the  high  order  control  system  components,  to  the  pitch  acceleration  at  time  t  =  04  excluding  these 
same  control  components,  following  a  step  control  input. 

Comparing  equations  (7)  and  (8)  it  can  be  seen  that  the  equivalent  system  CAP  will  not,  in 
general,  be  the  same  as  the  high  order  system  CAP'.  Although  the  equivalent  system  parameters 
can  be  correlated  unto  themselves,  they  will  not  correlate  with  the  high  order  system  parameters 
which  they  are  being  touted  to  represent.  However,  it  was  noted  in  reference  (q)  that  the  equivalent 
pitch  acceleration  has  characteristics  similar  to  Di  Franco's  nondimensionalizing  pitch  acceleration. 
Therefore,  an  equivalent  system  attenuation  factor  can  be  defined  as  the  ratio  of  the  maximum 
pitch  acceleration  of  the  high  order  system  to  the  initial  acceleration  of  the  low  order  system,  which 
will  correlate  the  two  representations.  It  was  also  shown  in  reference  (s),  using  available  inflight 
research  results,  that  the  CAP'  boundaries  of  table  X  can  be  utilized  to  correlate  high  order  re¬ 
sponses  with  pilot  opinion  ratings. 
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TABLE  X 

ATTENUATED  CONTROL  ANTICIPATION  PARAMETER  BOUNDARIES  -  CATEGORY  A 


Level 

CAP' 

Min 

Max 

1 

0.25 

1.5 

2 

0.15 

— 

Two  examples  will  serve  to  illustrate  the  differences  in  utilizing  CAP'  as  opposed  to  CAP  to 
correlate  the  equivalent  short  period  responses.  The  parameters  of  interest  for  two  selected  cases 
are  summarized  in  table  XI. 


TABLE  XI 

CAP  AND  CAP'  COMPARISONS 


Parameter 

A-6  Airplane  —  0/5eD 
0.72M/20.000  ft  altitude 

F-14  Airplane  —  0/F 
0.5M/15,000  ft  altitude 

LOS 

LOS 

HOS 

Ly 

“FIXED 

l“free 

HOS 

L„ 

“FIXED 

l“free 

“sp 

4.86 

4.75 

4.84 

2.78 

1.74 

2.88 

n/a 

24.8 

25.0 

31.2 

12.6 

12.7 

73.6 

e  (o+) 

0 

0.507 

0.444 

0 

0.0278 

0.0172 

^max 

0.305 

0.507 

0.444 

0.021 

0.0278 

0.0172 

D 

N 

CO 

CO 

0.612 

0.561 

0.590 

0.141 

0.117 

0.153 

wsp  2/n/“ 

0.952 

0.903 

0.751 

0.613 

0.238 

0.113 

CAP  =  0(O+)/nZss 

0 

0.904 

0.753 

0 

0.238 

0.113 

CAP' 

0.498 

0.544 

0.517 

0.150 

0.182 

0.139 

From  the  high  order  system  data,  it  can  be  seen  that  the  classical  definition  of  CAP  provides 
no  information  about  the  resulting  response  since  the  initial  acceleration  is  zero.  If  the  aircraft 
short  period  root  is  utilized  to  calculate  CAP,  false  indications  again  result  since  the  attenuation 
effects  of  the  control  system  are  not  accounted  for.  Additionally,  the  equivalent  system  models 
yield  differing  CAP  values  for  the  La-fixed  and  free  cases,  neither  of  which  are  equal  to  o?Sp^/n/o: 
of  the  high  order  system. 

The  premise  under  which  the  short  period  frequency  requirements  are  established  indicate 
that  as  acceleration  sensitivity  is  increased,  short  period  frequency  should  also  be  increased  to  main¬ 
tain  constant  pilot  rating  (i.e.  constant  control  anticipation  parameter).  It  could  therefore  be  as¬ 
sumed  that  for  the  low  order  system  to  be  rated  by  a  pilot  as  being  equivalent  to  the  high  order 
system,  this  same  relationship  should  hold.  However,  as  La  is  freed  in  the  search  routine,  the  CAP 
value  is  reduced.  Such  observations  would  tend  to  discount  the  equivalent  system  procedures. 
However,  if  the  differences  in  the  pitch  acceleration  characteristics  are  accounted  for  via  an  attenu¬ 
ation  factor,  very  good  agreement  is  found  among  the  example  cases. 
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Referring  back  to  figures  7  and  8,  there  is  no  constant  relationship  evident  between  each  of  the 
system  models  and  CAP.  The  data  necessary  to  define  CAP  and  CAP'  for  each  of  the  conditions 
analyzed  in  this  report  are  presented  in  tables  XI!  and  XIII.  Very  good  agreement  between  CAP' 
and  CAP'e  is  seen  to  exist  for  all  cases  with  the  largest  variation  in  the  two  parameters  arising  from 
the  configurations  with  the  highest  mismatch. 

The  data  may  be  viewed  in  the  format  of  Ml  L-F-8785C  by  determining  an  attenuated  or  effec¬ 
tive  frequency  for  which  CAP'  is  constant.  From  equation  (8)  this  effective  frequency  can  be  de¬ 
fined  as: 


^EFF  =  wsp  VTT  =  V(CAp,)  (n/a)  (9) 

This  parameter  is  most  easily  obtained  for  the  high  order  system  from  the  relationship: 


while  for  the  equivalent  systems 


/0maxHOS 

-EFF  *  -spe 

This  effective  or  attenuated  short  period  frequency  is  plotted  as  a  function  of  n/a  in  figures  9  and 
10.  The  level  1  and  2  boundaries  for  configuration  CR  are  those  presented  in  table  X.  CAP'  bound¬ 
aries  have  not  been  determined  for  approach  configurations.  Reference  (d)  indicates  that  Category 
A  (maneuvering)  boundaries  are  applicable  to  the  touchdown  portion  of  a  landing,  and  Category  C 
(approach)  boundaries  are  applicable  to  the  approach  phase.  Therefore,  as  a  point  of  reference  in 
defining  constant  CAP',  the  Category  A  boundaries  of  table  X  were  included  in  figure  10. 

The  data  presented  in  figures  9  and  10  not  only  indicate  that  CAP'  provides  an  additional 
measure  of  consistency  of  equivalent  systems,  but  that  CAP'  provides  a  more  accurate  representa¬ 
tion  of  frequency  characteristics  for  advanced  aircraft  than  does  CAP  The  effective  frequency  data 
for  each  flight  condition  analyzed  lies  along  a  line  of  nearly  constant  CAP'  for  all  of  the  various 
search  methods  (including  the  dominant  root  analysis).  As  a  result,  a  particular  level  of  flying 
qualities  can  now  be  identified  for  each  flight  condition.  (Unfortunately,  no  pilot  opinion  data 
obtained  during  maneuvering  tasks  is  available  to  verify  the  identified  levels). 

Comparing  figures  7  and  9  and  figures  8  and  10,  interesting  trends  can  be  noted  in  the  identi¬ 
fied  level  of  flying  qualities.  Those  aircraft  for  which  the  control  system  does  not  significantly 
affect  the  frequency  characteristics  (A-6  and  S-3)  indicate  solid  level  1  flying  qualities  whether  com¬ 
pared  against  CAP  or  CAP'.  However,  the  other  aircraft,  in  which  the  control  system  attenuates  the 
response,  are  nearer  the  lower  frequency  boundary  when  correlated  with  CAP'  rather  than  CAP. 
Since  CAP'  is  calculated  from  the  aircraft  responses  (pitch  and  normal  acceleration),  it  is  inter¬ 
preted  as  more  accurately  representing  the  aircraft's  characteristics. 

The  utilization  of  CAP'  also  provides  the  opportunity  for  analyzing  the  results  of  other  than 
first  over  second  order  equivalent  systems  with  regard  to  specification  requirements.  It  was  shown 
in  reference  (s)  that  for  systems  with  higher  order  poles  and  zeros,  CAP'  can  be  defined  as: 
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Notes:  (1)  Phugoid  contribution  included  in  high  order  system 

(2)  Configuration  PA(1) 

(3)  Configuration  PA(2) 


ATTENUATED  UNDAMPED  NATURAL  FREQUENCE 


NADC-81 069-60 


A-6  AIRPLANE 

O  .AN 

29,000  FT  ALTITUDE  □  .72H 
A  .88N 


S-3  AIRPLANE 

E  5,00ft  FT  ALTITUDE  ® 

li  .jin 


P-18  AIRPLANE 

10,000  FT  ALTITUDE  O  .SN 


Figure  9  —  CAP'  Correlation  of  Cruise  Configuration  Short  Period  Frequency  Characteristics 
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Figure  9  (continued)  —  CAP'  Correlation  of  Cruise  Configuration  Short  Period  Frequency 
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Frequency  Characteristics 
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CAP'e 


^sp2  gmax  H0S  l  pe 
n/a  9  (t  =  r)L0S’7rze 


(12) 


where  Pe  and  ze  represent  the  poles  and  zeros,  respectively,  necessary  to  improve  the  response 
match  in  the  frequency  range  of  interest.  Equation  (12)  was  utilized  to  compute  CAP'e  for  the 
F-14  and  A-7  cases  in  which  first  over  third  order  equivalent  systems  were  determined.  The  results 
for  each  of  these  cases  are  included  in  tables  XI I  and  XI 1 1.  CAP'e  for  each  of  the  first  over  third 
order  matches  showed  better  correlation  with  the  high  order  system  CAP'  than  did  the  first  over 
second  order  equivalent  systems.  This  improvement  in  CAP'e  can  also  be  correlated  directly  with 
the  mismatch  parameter  —  as  mismatch  is  improved  so  is  the  high  order  vs  low  order  CAP'.  The 
attenuated  control  anticipation  parameter  can,  then,  be  used  as  a  method  of  correlating  increased 
order  equivalent  systems  with  specification  frequency  requirements. 


CONCLUSIONS 


Longitudinal  short  period  equivalent  models  have  been  determined  which  represent  the  high 
order  dynamics  of  five  tactical  Navy  aircraft.  The  frequency  response  matching  method  utilized 
was  straightforward  and  easy  to  implement.  The  methods  applied  in  determining  the  best  low  order 
system  matches  can  be  divided  into  three  categories  based  on  their  frequency  content:  low, 
medium,  and  high. 

•  Low  frequency  components  were  essentially  ignored  in  the  matching  process  Where  the 
identified  equivalent  short  period  roots  were  closely  coupled  to  the  ignored  low  fre¬ 
quency  roots  (power  approach  configurations)  it  was  necessary  to  completely  remove 
their  contribution  from  the  high  order  representation  in  order  to  obtain  the  best  modal 
parameter  correlation. 

•  High  frequency  components  were  adequately  modelled  by  a  time  delay  parameter. 

•  Mid  frequency  components,  in  the  pilot's  frequency  range  of  interest,  compound  the 
matching  process.  When  numerator  and  denominator  roots  of  approximately  equal 
magnitude  are  introduced  in  the  high  order  system,  the  matching  process  easily  identifies 
a  set  of  equivalent  system  roots.  When  separation  occurs  between  the  numerator  and  de¬ 
nominator  roots,  or  an  unequal  number  of  roots  are  included,  the  frequency  matching 
technique  identifies  an  apparent  numerator  root  for  the  classical  first  over  second  order 
short  period  approximation,  which  may  be  significantly  different  from  the  of  the  air¬ 
plane.  This  condition  may  be  alleviated  by  either  1)  simultaneously  matching  pitch  rate 
and  normal  acceleration  and  accepting  the  (relatively)  large  mismatch  or  2)  introducing 
additional  roots  into  the  equivalent  system  model  to  account  for  the  dominating  high 
order  system  components. 

Correlation  of  the  equivalent  system  modal  parameters  with  Ml  L-SPEC  requirements  was 
straightforward  for  both  short  period  damping  ratio  and  time  delay.  Correlation  of  short  period 
frequency  characteristics  for  those  aircraft  with  significant  control  compensation,  however,  showed 
variation  in  identified  flying  qualities  levels.  The  level  of  flying  qualities  for  the  high  order  system 
was  different  from  the  equivalent  solution  with  L^-fixed  and  from  the  La-free  solution  when  anal¬ 
yzed  via  the  traditional  control  anticipation  parameter.  These  results  led  to  the  definition  of  a  con¬ 
trol  anticipation  parameter  which  allows  the  correlation  of  both  high  order  and  equivalent  low 
order  systems  to  pilot  opinion  ratings  and  hence  to  specification  boundaries.  This  attenuated  con¬ 
trol  anticipation  parameter  not  only  correlates  the  classical  short  period  equivalent  system  model 
with  the  physical  characteristics  of  the  high  order  system,  but  provides  the  opportunity  for  com¬ 
paring  alternate  low  order  equivalent  system  forms  with  flying  qualities  requirements. 


36 


NADC-81069-60 


The  A-6,  A-7,  S-3,  and  F-18  aircraft  responses  all  resulted  in  level  1  flying  qualities  for  each  of 
the  parameters  analyzed,  at  all  flight  conditions.  The  A-7  frequency  response  was,  however,  only 
marginally  level  1,  with  a  number  of  data  points  lying  directly  on  the  boundary.  The  F-14  airplane 
exhibited  level  1  performance  as  indicated  by  the  response  to  control  position  commands.  However, 
in  response  to  force  inputs,  the  F-14  airplane  would  be  predicted  to  be  a  level  2  airplane  both  in 
cruise  and  non-DLC  power  approach  conditions.  This  characteristic  is  evidenced  by  the  high  time 
delay  obtained  from  the  first  over  second  order  equivalent  analysis  and  from  the  effective  frequency 
characteristics  as  correlated  with  CAP'.  This  apparent  degradation  in  the  F-14's  handling  qualities 
can  be  traced  to  the  feel  system  implementation. 

RECOMMENDATIONS 

Frequency  response  matching  techniques  provide  an  effective  means  of  obtaining  low  order 
equivalent  models  of  high  order  longitudinal  aircraft  responses.  When  using  the  classical  short 
period  approximation  as  the  equivalent  system  model,  pitch  rate  and  normal  acceleration  should 
be  matched  simultaneously,  constraining  the  characteristic  denominator  roots  to  be  equal,  to  reduce 
the  variation  in  La  —  the  numerator  root. 

The  ratio  of  the  maximum  pitch  acceleration  (rather  than  the  initial  short  period  component) 
to  the  steady  state  normal  acceleration  (CAP')  experienced  in  a  pullup  should  be  utilized  as  the  cor¬ 
relating  parameter  for  Ml  L-SPEC  frequency  characteristics.  Additional  work  should  be  performed 
to  verify  the  Category  A  CAP'  boundaries  utilized  in  this  report  and  to  establish  similar  Category  C 
boundaries. 
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Data  describing  each  of  the  airplanes  and  their  respective  control  systems  were  obtained  from 
available  aerodynamic  stability  and  control  reports,  references  (h)  through  (m).  This  appendix  briefly 
describes  the  subject  airplanes  and  presents  a  block  diagram  of  their  respective  longitudinal  control 
system  as  modelled  in  this  analysis. 

A-6  —  The  A-6  airplane  is  a  twin  turbojet,  land  and  carrier  based,  subsonic,  all-weather  attack  air¬ 
craft.  Longitudinal  control  is  transmitted  from  the  pilot's  control  stick,  via  bellcranks  and  pushrods, 
to  an  all-moving  horizontal  stabilizer.  The  control  stick  feel  system  and  bobweight  arrangement  was 
not  included  in  the  present  model.  The  basic  aircraft  stability  is  augmented  by  the  feedback  of 
washed  out  pitch  rate  to  the  horizontal  stabilizer.  A  simplified  block  diagram  of  the  A-6  airplane's 
longitudinal  control  system,  as  modelled  in  this  analysis  is  presented  in  figure  A-1. 


A-6 


Figure  A-1  —  A-6  Airplane  Longitudinal  Control  System  Block  Diagram 

The  A-6  airplane’s  pitch  rate  response  to  pilot  longitudinal  control  inputs  can  be  represented 
by  a  fourth  order  numerator  over  sixth  order  denominator  transfer  function  of  the  form: 

» (si _ »n>  N<?  i2s+  ^  (a_, 

Sep(s)  A  (2s  +  111.03s  +  1)  +  N  j  (.24s) 

where  JN  =  0.0122  rad/in  in  cruise  configuration  and  0.0232  rad/in  in  power  approach  configura¬ 
tion.  The  transfer  functions  representing  the  A-6  airplane,  as  analyzed  in  this  report,  are  presented 
in  table  A-!. 

S-3  —  The  S-3A  airplane  is  a  twin  turbofan  powered,  land  and  carrier  based,  subsonic,  anti¬ 
submarine  warfare  aircraft.  Longitudinal  control  is  accomplished  via  a  mechanical  control  system 
which  operates  the  elevator.  Control  stick  dynamics  were  not  included  in  the  model  of  the  S-3 
aircraft.  There  is  no  longitudinal  stability  augmentation  system  included  in  the  aircraft.  A  block 
diagram  of  the  S-3's  control  system,  as  modelled  in  this  analysis,  is  presented  in  figure  A-2. 
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Figure  A-2  —  S-3  Airplane  Longitudinal  Control  System  Block  Diagram 


A-2 


NADC-81 069-60 


TABLE  A-l.  A-6  AIRPLANE  TRANSFER  FUNCTIONS 


Config¬ 

uration 

Altitude 

(ft) 

Airspeed 

(M/KEAS) 

Transfer  Function 

CR 

20,000 

0.4/179 

6 

4.31  (0)  (.0147)  (.506)  (.5) 

5ep 

[0.029,  .11]  [0.63,2.32]  (0.499)  (31.96) 

CR 

20,000 

0.72/323 

9 

13.94  (0)  (0.011)  (1.077)  (0.5) 

5ep 

[0.088,0.043]  [0.86,4.86]  (0.428)  (28.12) 

CR 

20,000 

0.88/395 

d 

19.99  (0)  (0.028)  (1.341)  (0.5) 

5ep 

(0.069)  (-0.045)  [0.93,6.82]  (0.415)  (24.68) 

PA 

Sea  Level 

0.14/95 

B 

4.26  (0)  (0.186)  (0.461)  (0.5) 

Sep 

[0.048,0.26]  [0.71,1.46]  (0.5)  (32.63) 

0<1>  

4.26  (0.578)  (0.5) 

Sep 

[0.70,  1.47]  (0.481)  (32.63) 

Note:  (1)  Phugoid  contributions  ignored. 


The  S-3  airplane's  pitch  rate  to  cockpit  control  position  command  may  be  represented  by  a 
third  order  numerator  over  fifth  order  denominator  transfer  function  of  the  form: 

9Js)  m  (0-0052)  Ng(s)  1  (a_2) 

5ep(s)  Ms)  '  0.0294s  +  1 

The  transfer  functions  representing  the  S-3  airplane,  as  analyzed  in  this  report,  are  presented  in 
Table  A-ll. 

F-14  —  The  F-14  airplane  is  a  twin  turbo-fan  powered,  land  and  carrier  based,  supersonic  fighter 
aircraft.  Longitudinal  control  is  accomplished  via  an  irreversible  mechanical  flight  control  system 
which  transmits  cockpit  control  commands  to  an  all  moving  horizontal  stabilizer.  Fore  and  aft 
bobweights,  not  modelled  in  this  analysis,  are  utilized  to  provide  tailored  stick  force  per  g  character¬ 
istics.  Pilot  force  feel  is  provided  via  a  nonlinear  feel  spring  and  stick  motions  are  damped  by  a 
sprashpot.  The  airplane's  basic  stability  is  augmented  through  the  feedback  of  washed  out  pitch 
rate  which  is  fed  through  a  shaping  network  to  obtain  the  desired  response.  A  block  diagram  of  the 
F-14's  longitudinal  control  system,  as  modelled  in  this  analysis,  is  presented  in  figure  A-3.  The  F-14 
airplane's  response  to  longitudinal  control  inputs  may  be  represented  by  the  following  general 
transfer  function: 

X(s)  (0.0374)  Nx  (2s  +  1 )  (0.53s  +  1 )  (0.071 5s  +  1 ) 

-r -  - - - -  (a-3) 

6ep(s)  A(0.05s  +  1)  (2s  +  1)  (0.53s  +  1)  (0.0715s+  1)  +2Kqs  (0.2s  +  I)2  N^ 

The  transfer  functions  representing  the  F-14's  response  to  cockpit  control  position  inputs,  as 
analyzed  in  this  report,  are  presented  in  table  A-l  II. 
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TABLE  A-l I 

S-3  AIRPLANE  TRANSFER  FUNCTIONS 


Config¬ 

uration 

Altitude 

(ft) 

Airspeed 

(M/KEAS) 

Transfer  Function 

CR 

15,000 

0.36/179 

6 

249.2  (0)  (0.0227)  (0.714) 

5eP 

[0.048,  .106]  [0.44,2.59]  (34.01) 

CR 

1 5,000 

0.71/353 

6 

786.7  (0)  (0.032)  (1.766) 

5eP 

[0.8,  .019]  [0.48,5.45]  (34.01) 

PA 

Sea  Level 

0.15/97 

6  _ 

73.43  (0)  (0.173)  (0.594) 

5eP 

[0.18,0.197]  [0.57,1.53]  (34.01) 

g(D_ 

73.43  (0.675) 

5ep 

[0.57,1.52]  (34.01) 

Note:  (1)  Phugoid  contributions  ignored. 


TABLE  A-l II 


F-14  AIRPLANE  TRANSFER  FUNCTIONS 


Config¬ 

uration 

Altitude 

(ft) 

Airspeed 

(M/KEAS) 

Transfer  Function 

CR 

1 5,000 

0.5/248 

6 

5.26  (0)  (0.0103)  (0.773)  (0.5)  (1.887)  (13.986) 

5eP 

[.016,  .082]  [.61,2.78]  (.418)  (1.34)  [.97,17.04] 

Nzcg  

-1.50  (0)  (0.00066)  (6.619)  (-6  728)  (0.5)  (1.887)  (13.986) 

5ep 

dC) 

NzCr  _ 

1.34  (0)  (0.00066)  (49.99)  (0.5)  (1.887)  (13.986) 

SeP 

D 

CR 

1 5,000 

0.7/348 

6 

11.48  (0)  (0.0119)  (0.962)  (0.5)  (1.887)  (13.986) 

5eP 

[0.057,  .065]  [0.58,4.27]  (0.408)  (1.532)  [0.91,17.51] 

Nzcg  

-3.22  (0)  (0.0079)  (8.816)  (-8.919)  (0.5)  (1.887)  (13.986) 

Sep 

D 

Nzcr  _ 

3.45  (0)  (0.0079)  (73.55)  (0.5)  (1.887)  (13.986) 

Sep 

D 

Note:  (1)  D  =  Denominator  of  0/5ep  transfer  function. 
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TABLE  A  ll  I 

F-14  AIRPLANE  TRANSFER  FUNCTIONS 
(Continued) 


Config¬ 

uration 

Altitude 

(ft) 

Airspeed 

(M/KEAS) 

Transfer  Function 

CR 

1 5,000 

0.83/412 

d 

15.22  (0)  (0.012)  (0.965)  (0.5)  (1.887)  (13.986) 

Sep 

[0.11,0.058]  [0.62,5.01]  (0.418)  (1.562)  [0.87,17.95] 

Nzcg 

-4.16  (0)  (0.0093)  (9.78)  (-9.856)  (0.5)  (1.887)  (13.986) 

S  " 

D 

Nzcr 

4.35  (0)  (0.0093)  (92.38)  (0.5)  (1.887)  (13  986) 

S 

D 

CR 

1 5,000 

1 .2/596 

d 

30.46  (0)  (0.0295)  (1.143)  (0.5)  (1.887)  (13.986) 

Sep 

[0.25,0.064]  [0.71,6.88]  (0.418)  (1.685)  [0.70,20.37] 

Nzcg 

-7.04(0)  (0.0278)  (13.75)  (-14.17)  (0.5)  (1.887)  (13.986) 

Sep  ’ 

D 

Nzcr_ 

11.02  (0)  (0.0278)  (124.4)  (0.5)  (1.887)  (13.986) 

Sep 

D 

PA(1)(2) 

Sea 

Level 

.19/126 

Q 

1.179  (0.473)  (0.5)  (1.887)  (13.986) 

Sep  ' 

[0.52,  1.48]  (0.506)  (1.591)  (15.09)  (18.66) 

Nzcg  

-5.587  (3.185)  (-2.735)  (0.5)  (1.887)  (13.986) 

Sep 

D  (D 

NzCr 

0.134(363.6)  (0.5)  (1.887)  (13.986) 

Sep 

D 

PA(2)(2) 

Sea 

Level 

.18/121 

6 

1.034  (0.444)  (0.5)  (1.887)  (13.986) 

Sep 

(0.7,  1.05)  (0.531)  (1.48)  (14.91)  (18.87) 

Nzcg_ 

-5.899  (-2.357)  (2.756)  (0.5)  (1.887)  (13.986) 

5eP 

D 

NzCr  

0.264  (145.49)  (0.5)  (1.887)  (13.986) 

Sep 

D 

Notes:  (1)  D  =  Denominator  of  0/5 6p  transfer  function. 
(2)  Phugoid  contributions  ignored. 
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Figure  A-3  —  F-14  Airplane  Longitudinal  Control  System  Block  Diagram 
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The  response  to  cockpit  force  commands  is  obtained  from: 


X(s)  =  X(s)  5eP(s) 
Fx(s)  5ep(s)  Fx(s) 


(a-4) 


where: 


S(s) 

Fx(s) 


- 1 - - - 1 -  ,  5en<  1  'n 

(s  +  3.366)  (sz  +  36.45s  +  1580)  p 


26.825  (s  +  39.815) 


(a-5) 


A-7  —  The  A-7  airplane  is  a  single  place  turbo  fan  powered,  land  and  carrier  based,  light  attack  air¬ 
craft.  It  contains  an  irreversible  mechanical  longitudinal  control  system  with  both  stability  and  con¬ 
trol  augmentation.  The  stability  augmentation  system  provides  pitch  rate  and  filtered  normal  ac¬ 
celeration  feedback  signals  to  augment  the  aircraft's  basic  stability  characteristics.  The  command 
augmentation  system  feeds  control  force  signals  forward  through  a  prefilter  as  a  means  of  increasing 
the  pilot's  commanded  input.  The  A-7  control  stick  dynamics  were  not  modelled  in  this  analysis. 

A  simplified  block  diagram  of  the  A-7  airplane's  longitudinal  control  system  is  presented  in  figure 
A-4.  The  A-7  airplane's  response  to  control  force  inputs  may  be  represented  by  the  following  gen¬ 
eral  transfer  function: 


NX  [(0.55s  +1)Km  +  Kcas] 


X(s) 


(a-6) 


Fx(s) 


A  (0.55s  +  1 )  (0.05s  +  1 )  +  KnZa  NnZa  +  * 


The  transfer  functions  representing  the  A-7  airplanes  response  to  pilot  force  commands,  as  modelled 
in  this  analysis,  are  presented  in  table  A-IV. 

F-18  —  The  F-18  is  a  single-place,  turbo-powered,  land  and  carrier  based,  fighter  aircraft  controlled 
by  a  digital  flight  control  system.  Separate  flight  control  laws  are  provided  for  differing  flight 
regimes.  Electrical  signals  are  generated  from  the  pilot's  control  force  inputs,  passed  to  the  com¬ 
puter,  modified  by  various  gain  and  shaping  networks,  and  finally  passed  to  the  all-moving  hori¬ 
zontal  stabilizer.  Normal  acceleration  and  pitch  rate  signals  are  also  input  to  the  computer,  where 
they  are  shaped  and  gain  scheduled  before  being  summed  with  the  command  input  signals.  The 
cockpit  control  feed  system  dynamics  were  not  included  in  this  analysis.  Simplified  block  diagrams 
of  the  F-18  airplane's  longitudinal  flight  control  system  (revision  4.1)  are  presented  in  figures  A-5 
and  A-6  for  the  Cruise  and  Power  Approach  flight  regimes,  respectively.  The  F-18  high  order  sys¬ 
tem  transfer  functions  investigated  in  this  analysis,  obtained  from  the  NAVAI  RDEVCEN  Flight 
Control  Section,  Code  6012,  are  presented  in  table  A-V. 
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Figure  A-4  —  A-7  Airplane  Longitudinal  Control  System  Block  Diagram 
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TABLE  A-l V 

A-7  AIRPLANE  TRANSFER  FUNCTIONS 


Config¬ 

uration 

Altitude 

(ft) 

Airspeed 

(M/KEAS) 

T ransfer  Function 

CR 

15,000 

0.3/149 

6 

0.158  (0)  (-0.0082)  (0.506)  (7.272) 

F 

[0.064,0.12]  [0.52,1.83]  (2.079)  (18.758) 

n7 

zcg  

-0.833  (-0.0042)  (-0.05)  (5.55)  (-5.08)  (7.272) 

F 

D^ 

n7 

zcr  

0.08  (-0.0042)  (-0.05)  (340)  (7.272) 

F 

D 

CR 

15,000 

0.6/298 

6 

0.658  (0)  (0.0072)  (1.09)  (7.272) 

F 

[0.068,  .05]  [0.76,3.23]  (7.069)  (1 1.84) 

CR 

15,000 

0.9/447 

6 

1.456  (0(  (0.0443)  (1.97)  (7.272) 

F 

(0.0644)  (-0.0255)  [0.78,4.01]  [0.64,  14.94] 

Note:  (1)  D  =  Denominator  of  0/56p  Transfer  Function 


TABLE  A-V 

F-18  AIRPLANE  TRANSFER  FUNCTIONS 


Config¬ 

uration 

Altitude 

(ft) 

Airspeed 

(M/KEAS) 

Transfer  Function 

CR 

10,000 

0.5/274 

_s_ 

F 

6284.5(0X0.25)  (1X1 )  (1 .247)  (1 .25)  (3.33X30.03)  [0.03,60] 

(-0.00261(0.22)  (0.81, 1.11]  (0.89X1. 67)  (0.76, 4.47]  [0.76,30]  (0.35,55]  [-0.80,61.7] 

PA 

Sea 

Level 

.20/133 

& 

F 

170.7(0.395X1  )(1 1(14.131(20 .321(25.08)  [0.03,60] 

[0.682,1 .78]  (0.368)  (2X1 5.82)  (19.82X24.71 )  [0.88,50]  [0.66,65.1  ] 

nzCr 

F 

2.4 5(  1 X 1  X8.85)  [0.88,22.7]  (32.06)  [0.03,60] 

D(1) 

Note:  (1)  D  =  Denominator  of  0/6ep  Transfer  Function 
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Figure  A-5  —  F-18  Airplane  Longitudinal  Control  System  Block  Diagram  —  Cruise  Configuration 
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Figure  A-6  —  F-18  Airplane  Longitudinal  Control  System  Block  Diagram  —  Power  Approach  Configuration 


APPENDIX  B 


Frequency  and  Time  History  Response  Comparisons 


Bode  plots  showing  frequency  response  and  time  history  plots  showing  the  airplane's  calcu¬ 
lated  response  to  step  control  inputs  are  compared  for  the  assumed  high  order  and  various  low 
order  systems  in  the  illustrations  of  this  appendix,  figures  B-1  through  B-24.  Data  are  presented  for 
each  of  the  airplanes  analyzed  at  one  flight  condition  in  both  cruise  and  power  approach  configura¬ 
tions. 
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Figure  B-24  —  F-18  Airplane  Normal  Acceleration  Response  -  Configuration  PA  -  Simultaneous  Q  and  Nz  Matching 
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